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ABSTRACT 

We present the results of a study comparing the stellar populations in the elliptical galaxies of 
Hickson Compact Groups (HCGs) with those in low density environments. We analyze a sample of 
intermediate resolution (4.25A FWHM) spectra of 22 galaxies in HCGs and 12 in the field or in 
loose groups. Three different population synthesis models and stellar population analyses are used to 
make the results more robust. Great care has been taken to correct for the emission contamination 
of the Balmer lines, used to determine the mean luminosity weighted ages of the galaxies. We find 
that, on average, galaxies in HCGs are systematically older and more metal-poor than in the field, 
in agreement with previous works. The most interesting finding is that the low-cr galaxies in HCGs 
show an enhanced [Mg/Fe] ratio and depleted metallicity [Z/H] with respect to their counterparts in 
the field. This anomalous behavior is interpreted as evidence for the action of a mechanism which 
truncated the star formation (SF). Hydrodynamical simulations of galaxy mergers (Di Matteo et al. 
2005) support this interpretation by predicting the quenching of star formation soon after the merger 
event. Combining this scenario and the evidence presented here, the HCGs, generally considered to 
be ideal environments for galaxy-galaxy interactions, become ideal places for SF truncation. 
Subject headings: galaxies: Compact groups - galaxies: Evolution - galaxies: Interactions - galaxies: 
Clustering ~ galaxies: Stellar Populations 



1. INTRODUCTION 

In the last decade, work on Hickson Compact Groups 
(HCGs) has been characterized by a sense of disappoint- 
ment. On one hand, HCGs of galaxies, with their high 
spatial densities and low velocity dispersions appear to 
be ideal places for interactions and mergers. On the other 
hand, the interaction-activity paradigm (e.g. Barnes & 
Hernquist 1992) predicts the triggering of star formation 
and a widespread flux increase at almost every wave- 
length as a result of these interactions. Although ob- 
servations have found abundant traces of tidal interac- 
tions in HCG galaxies (Mendes de Oliveira & Hickson 
1994), their star formation levels are surprisingly simi- 
lar to those found in isolated galaxies (e.g. Zepf, Whit- 
more & Levison 1991; Moles et al. 1994; AUam et al. 
1996; Verdes-Montenegro et al. 1998; Iglesias-Paramo & 
Vi'lchez 1999). 

There is evidence that this unexpected behavior is not 
exclusive to the HCGs, but is shared by the broader fam- 
ily of galaxy groups and clusters. The exploitation of 
large scale redshift surveys has allowed the study of large 
samples of galaxies in clusters (Hashimoto et al. 1998) 
and groups (Wilman et al. 2005). One of the most ap- 
pealing results of these works is the realization that the 
dense galactic environments, found in groups and clus- 
ters, are hostile to star formation. Unfortunately, the 
mechanism which truncates the star formation in groups 
remains virtually unknown. Although, there is no short- 
age of proposed mechanisms to halt star formation in 
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dense environments, they are considered ineffective in 
group conditions. In addition, the merger of galaxies, fa- 
vored by the conditions in HCGs, has also been usually 
excluded because it enhances the star formation instead 
of stopping it. 

This view has been modified based on recent hydro- 
dynamical simulations of galaxy mergers by Di Matteo, 
Springel & Hernquist (2005), which take the feedback 
from supermassive black holes into account. They show 
that the star formation enhancement is just a transient 
phase of a longer process which, ends up depleting the 
gas and quenching both star formation and nuclear activ- 
ity. A short time (~ 0.5 Gyr) after the merger event, the 
remnant has evolved to a dead spheroid where the tracers 
of the transient starburst begin to fade. The short-lived 
SF tracers, such as the Ha emission, disappear along 
with the gas, while the color, which slowly reddens after 
star formation has ceased, reaches u-r ~ 2.2-2.3. This 
places the remnant in the red sequence of the bimodal 
color distribution in about 1-2 Gyr (Springel, Di Matteo 
& Hernquist 2005). From this new perspective, the fail- 
ure of HCGs to display traces of ongoing or recent star 
formation would not seem so disappointing. 

A more appropriate tool to determine the star for- 
mation histories (SFH) of galaxies is the study of their 
stellar populations. The use of this technique has pro- 
duced valuable insights into the history of large sam- 
ples of galaxies in high and low-density environments 
(Kuntschner et al. 2002; Thomas et al 2005). There 
also exist a few stellar population studies with smaller 
samples of early-type galaxies in HCGs (do la Rosa et 
al. 2001b; Proctor et al. 2004; Mendes de Oliveira et al. 
2005). Their common result is that early-type galaxies 
in HCGs are older and more metal-poor than their coun- 
terparts in the field. Opposite results have also been 



2 



reported (Ferreras et al. 2006). Our goal in this work 
is to use stellar populations to interpret the relics of the 
past episodes of galaxy mergers. The present study su- 
persedes a previous work (de la Rosa et al. 2001b), using 
the same sample but with a less sophisticated analysis. 

The study of stellar populations in galaxies has reached 
maturity. It has overcome several obstacles, like the 
agc-mctallicity degeneracy, which was partially lifted 
with the introduction of suitable Lick spectral indices 
(Worthey 1994). The non-solar element abundance ra- 
tios in the galaxies was taken into account with the inclu- 
sion of the abundance ratios effects in the stellar popula- 
tion models (Weiss, Peletier & Mattcucci 1995; Vazdckis 
et al. 1997; Tantalo, Chiosi & Bressan 1998; Trager 
et al. 2000a; Thomas, Maraston & Bender 2003, here- 
after TMB03; Barbuy et al. 2003). Despite these suc- 
cesshil achievements, stellar population studies are still 
not problem-free. Although novel approaches and better 
models appear at a good pace, the convergence of their 
results is not accompanying these improvements. Differ- 
ent stellar population approaches applied to the same 
galaxy do not generally agree on the absolute values 
of the derived parameters. Another drawback of stel- 
lar population studies arises from the relatively frequent 
emission contamination of the Balmer absorption indices, 
which are extensively used to determine galaxy ages. In 
the present study, the emission fill-in of the Balmer lines 
is a potential source of systematic errors, due to a possi- 
ble different incidence of AGNs in the two environments 
(Coziol, lovino & de Carvalho, 2000). If left uncorrected, 
the emission fill-in weakens the Balmer absorption lines 
and leads us to derive systematically larger ages and 
lower metallicities. Some authors have abandoned emis- 
sion correction of their galaxies, simply excluding the 
contaminated cases from their samples. 

In this work we use medium-resolution spectra of a 
homogeneous sample of early type galaxies belonging to 
both HCGs and the field or loose groups. Our goal is to 
study the effect of the environment on the stellar popu- 
lation parameters (age, metallicity and abundance ratio) 
of galaxies within the sample. For this reason, the ab- 
solute values of these parameters are much less relevant 
than their relative values. We have taken special care to 
double-check the influence of both the choice of a partic- 
ular stellar population analysis and the emission contam- 
ination on the results. We have carried out three parallel 
analysis of the same data using very different approaches. 
We first apply the Lick procedure with the low resolu- 
tion TMB03 models, including the effects of non-solar 
element abundance ratios. We then use the SSP models 
presented in Mendes de Oliveira et al. (2005, hereafter 
MCGB05), which also work at the Lick resolution and 
include non-solar abundance effects. Because the work 
by Mendes de Oliveira et al. (2005) contains two kinds 
of data relevant for this study, i.e. SSP models as well 
as galaxy parameters, we use the acronym MCGB05 to 
reference the models, and the full reference to cite the 
galaxy results, especially in subsection 5.2. Third is the 
Vazdckis (1999, hereafter Vaz99) approach, with a high 
resolution Stellar Population (SSP) model and a different 
stellar population analysis strategy, where the models are 
adapted to the resolution of the data, instead of the usual 
converse method. Our aim is not to compare the perfor- 
mances of the three approaches, but to gain insight into 



the systematic differences to assess the robustness of the 
results. The emission correction has also been calculated 
following two distinct strategies. One of them, pioneered 
by Gonzalez (1993) (hereafter G93), calculates the H/3 
emission contamination from the emission of the nearby 
[OIII]A5007 line. The second approach, introduced by 
Caldwell, Rose & Concannon (2003), estimates the 11/3 
fill-in through several steps, involving the determination 
of the Ha emission and the use of the ratios between the 
Balmer line strengths in nebular emission spectra (Os- 
terbrock 1989). 

This paper is organized as follows: Section 2 describes 
the observations and the basic reduction process. In Sec- 
tion 3, the emission correction is discussed. Section 4 is 
devoted to the extraction of the stellar population pa- 
rameters with three alternative methods, whose results 
and correlations are shown in Section 5. In Section 6 a 
physical scenario is proposed to account for the results, 
and a summary is presented in Section 7. 

2. OBSERVATIONS AND BASIC REDUCTIONS 

2.1. The Data Sample 

Our sample consists of 22 bright elliptical galaxies lo- 
cated in the cores of HCGs (Hickson 1982) and 12 bright 

bona fide E/SOs located in the field, very loose groups 
or cluster outskirts, hereafter named the field galaxies. 
Long-slit observations were obtained during a four night 
run (1994 January 9-12) at the KPNO 2.1 m telescope 
with the GoldCam CCD spectrometer. By using a sin- 
gle telescope and optical configuration we guarantee the 
homogeneity of the sample. A 600 line/mm grating was 
used to yield a ~ 1.25 A /pixel dispersion and a 4.25 
A FWHM spectral resolution over the wavelength range 
3500 to 7000 A. A more detailed description of the ob- 
servations has been published in de la Rosa, de Carvalho 
& Zepf (2001a, hereafter Paper I)). 

The main details of the data have been summarized 
in Table 1, where column 1 shows the list of observed 
objects. In general, the slits were not aligned along a 
preferred axis of the galaxies. Repeated observations of 
the same galaxy were combined except for those cases 
where slit inclinations differed by more than 20 degrees. 
In those particular cases where several observations of the 
same object exist, each different slit inclination is marked 
with a plus (+) and treated as an independent object. A 
positive side effect of this discrimination is to increase 
the number of observations to populate the statistical 
relations. However the results of this study (Section 5) 
only include one observation per object, therefore, in the 
case of multiple inclinations only the observation with 
the highest S/N has been selected. 

The S/N per A shown in column 2, have been calcu- 
lated for the Re///2 aperture and the 4795 A to 5465 A 
spectral interval (RED). Calculations take into account 
Poisson statistics, with the gain (2.8 e~/ADU), readout- 
noise (8.5 e~), number of combined spectra, and the sky 
value. The S/N ranges from 20 to 250, with a median 
value of 55. We have not excluded any galaxies in the 
sample based on their S/N. 

2.2. Basic Reduction and Aperture Extraction 

Data reduction was performed using the standard 
IRAF tasks, including bias subtraction, flat-fielding, cos- 
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mic rays and bad pixel removal, etc. The extinction cor- 
rection and flux calibration of the spectra were made 
with the spectrophotometric standard stars Feige 34 and 
Hiltner 600, which were observed in each observing ses- 
sion. 

Aperture extraction was performed in a coherent way 
to allow for comparison of all galaxies in the sample and 
avoid the use of aperture corrections to the spectral in- 
dices. Extracted apertures are fractions of the effective 
radius (Re//); therefore, we have worked out the Re// 
value along the slit during the observations. Both the 
data for the calculation and the results are presented in 
Table 1. 

Isophotal photometry data, Reff (column 3) and el- 
lipticity (column 5) were taken from Zepf & Whitmore 
(1993) and Trager et al. (2G0Gb). Missing data for three 
of our sample galaxies (HCG 46a, HCG 51a and HCG 
62a) were estimated through indirect methods, using the 
roughly linear relation between the de Vaucouleurs diam- 
eter (-D25) and Reff - Ellipticities for these three galaxies 
were taken at isophote 25 from the Lyon-Meudon Extra- 
galactic Database (LEDA) 

The slit inclination is presented in column 4 as the 
angle between the major axis of the ellipse and the slit. 
This value has been calculated from the slit orientation 
and the value of the major axis position angle (PA) from 
LEDA. Three galaxies '(HCG 59b, HCG 93a and NGC 
4552), with rather low ellipticity, lacked the PA informa- 
tion, which were directly measured from images in the 
NASA/Ipac Extragalactic Database (NED). The previ- 
ous data are used to calculate the i?e/ / value along each 
particular slit inclination. Column 6 shows the aperture 
values for i?e///8 (in arcsecs) projected along the slit. 
Although only the Ref / /2 and R^ff/S apertures are used 
here, smaller fractions (1/4 or 1/6) replace the Reff/8 
apertures when their sizes are of the order of 3 arcsecs, 
the average seeing of the observations. For instance, a 
Ref //6 aperture would be preferred over a Ref / /8 one 
when both fall inside the same blurred bin, just because 
it contributes with more flux. Values for those cases are 
given in parentheses. 

2.3. Dynamical Parameters 

The spectral lines were corrected for two main sources 
of broadening: the rotation and the internal velocity dis- 
persion of the galaxy. The rotational correction is carried 
out along the aperture extraction, where every pixel-size 
slice is corrected for the rotational Doppler shift, before 
being stacked in the extracted spectrum. The Doppler 
shifts were measured with the FXCOR/IRAF task, via 
cross-correlation with a rest frame stellar template SAO- 
079251 (Paper I). The central velocity dispersions of the 
galaxies (ctq) have also been measured with the procedure 
used in Paper I, but averaging the results of three differ- 
ent variants: (i) using the SAO-079251 stellar template 
in the full spectrum; (ii) using the best fitting SSP-Vaz99 
models in the RED interval (4795 A to 5465 A ) and (iii) 
using the best fitting SSP-Vaz99 models in the BLUE 
interval (3856 A to 4476 A ). Values of ctq for Reff/8 are 
presented in column 2 of Table 2. 

3. EMISSION CORRECTION 

Residual emission in the Balmer line indices produce 
artificially large ages in the affected galaxies and rep- 
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resents a potential source of systematic errors for the 
present study. Extracting emission lines from an under- 
lying absorption spectrum is especially difficult for the 
Balmer lines, because they coincide with the absorption 
lines which are our main source of information about 
population ages. Due to this difficulty, we have double- 
checked the residual emission correction by using two 
alternative approaches. The first strategy, proposed by 
G93, uses metal forbidden lines, generally independent of 
any underlying absorption feature, to deduce the Balmer 
emission components. The second strategy, proposed 
by Caldwell, Rose & Concannon (2003), calculates the 
Ha emission to deduce the higher order Balmer emission 
lines. 



3.1. 



First Approach 

The first approach seeks a relation between two spec- 
tral indices: the UPemissiom with the same bandpasses 
defined for the standard UPabsorpUon, and the nearby 
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Fig. 1. — The use of two independent methods for correction of 
the up emission. In the first method (a) the relation between the 
emission in [OIII]A5007A and H/3 is used to derive the contamina- 
tion of the Balmer index. The second method (b) uses the relation 
between the Rose index Hn/Fe with Haabsorption to isolate the 
emission component of the Ha line, which is used to deduce the 
emission contamination of the H/3 line. Despite these rather scat- 
tered relations, the corrected H/3 indices, with the two methods, 
show a good match (c) . The dotted diagonal line in (c) corresponds 
to the equal-value line. 



[OIIIJA5007A index, whose bandpasses are defined in 
G93. Both indices are measured in the residual spec- 
tra, which are constructed by subtracting the best fitting 
SSP-Vaz99 model from the original spectrum. Only data 
showing emission in any of the indices were used for the 
correlation, resulting in a linear fit of UPemission = 0.180 
X [OIII]A5007 - 0.092, with rms=0.208 (Figure 1(a)). 



According to this statistical correction, an [OIIIJA5007 
emission of, say -0.4, implies a H/3 emission filling of - 
0.164, i.e. a AHP = 0.164 A correction to be added to 
the contaminated H/3 index. 

3.2. Second Approach and Comparison 

CaldwcJl. Rose & Coneannon's (2003) approach is ex- 
tensively explained in their paper and we will just outline 
the steps we have followed. In short, one measures the 
Ha emission to deduce the corresponding H/3 contam- 
ination. We divide our sample into two sets based on 
the presence or absence of [OIIIJA5007 emission. The 
template set is comprised of galaxies with a S/N greater 
than 40 and [OIII]A5007 > - 0.3, while the contaminated 
set contains galaxies with [OIII]A5007 emission (EW < 
- 0.3). We have used all our R/2 and R/8 spectra as 
independent data, resulting in 19 contaminated, and 48 
template spectra. Furthermore, to improve the quality 
of the template set, we have excluded seven spectra with 
the largest departures from the mean relation. Two addi- 
tional spectra were excluded due to poor spectral quality 
in the Ha region. We then measure the Ha absorption 
and the Rose index Hn/Fe for each template galaxy. A 
linear relation is found between Haabsorption and Hn/Fe 
(Figure 1(b)), suggesting that this last index is a reli- 
able predictor of the Ha absorption. Our least-squares 
fit gives Haabsorption = - 2.08 X (Hn/Fe) + 3.97, with 
a rms scatter of ± 0.19, in good agreement with the 
relation found by Caldwell, Rose & Concannon (2003). 
We make the assumption that this same linear relation 
also applies to the emission-contaminated galaxies, yield- 
ing an indirect method to estimate the HaabsorpUon by 
means of the Hn/Fe index. Each contaminated galaxy 
is assigned a template with both a similar Hn/Fe index 
and a lower velocity dispersion. The previous assump- 
tion means that, after the template has been broadened 
to match the a of the contaminated galaxy, both spec- 
tra should have similar Ha absorption. The Ha emission 
can now be measured in the residual spectrum (contam- 
inated - cr-broadened template). Once the Ha emission 
has been obtained and assuming the Balmer emission 
decrement ratios proposed by Osterbrock (1989), e.g. 
Ha/H/3 ^ 3.05, the higher order Balmer emission lines 
can be estimated. Note that the low dust extinction in 
our sample galaxies alleviates the need for corrections to 
the theoretical decrement ratios, because they are irrele- 
vant compared to the uncertainties of the models and the 
extinction curves. The intensity of emission in H/3 is not 
yet the correction we are seeking, because Lick indices 
do not only measure the spectral lines after which they 
have been named, but the neighboring lines included in 
the bandpasses. An artificial emission spectrum is cre- 
ated with the IRAF/mkldspec routine, including the Ha 
emission line, with the flux measured in the residual spec- 
trum, and the decremented higher order Balmer emis- 
sion lines. Finally, each artificial emission spectrum is 
subtracted from the contaminated one and the emission- 
corrected Balmer indices can be measured. 

Despite the fact that both methods are completely dif- 
ferent and based on rather scattered relations (Figures 
1 (a) and(b)), there is a good match between their re- 
sults. Figure 1(c) compares the emission-corrected H/3 
indices measured with the two alternative approaches. 
The comparison shows a rms scatter of 8%, with 4% 
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Fig. 2. — Comparison of the raw values for the most relevant 
spectral indices with G93. The dotted diagonal represents the 
equal-value line, while the average errors are plotted in the lower 
right corner. 
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systematically larger values of the corrected index ob- 
tained with the second approach. Translated into ages, 
the uncertainty and the shift correspond to roughly 2 
and 1 Gyrs, using the grid of Figure 5(c). Consider- 
ing the relatively good agreement between both correc- 
tions we have adopted their average as the final value. 
Only 19 contaminated spectra {Reff/8 and Reff/2) with 
EW[OIII]A5007 < 0.3, have boon corrected for emission 
and arc identified with an asterisk in Table 2. 

4. THE STELLAR POPULATION STUDY 

There are a variety of models and strategies to cal- 
culate stellar population parameters in galaxies. We use 
three alternative stellar population analyses with the aim 
of avoiding systematic effects due to the choice of models, 
rather than to evaluate the performance of each analysis. 
We first present the standard Lick/IDS approach, with 
the SSP models developed by TMB03, which include 
abundance ratio effects. A second approach, MCGB05, 
uses SSP models based on synthetic stellar spectra com- 
puted in Barbuy et al. (2003). The third approach is 
the model/analysis proposed by Vaz99, a novel strategy 
which leaves the original spectral resolution of the data 
untouched. 

4.1. The TMB03 Approach 

4.1.1. Transformation to Lick/IDS 

Although we lack standard stars in common with the 
Lick library, we have 12 galaxies in common with G93, 
one of the reference works in the Lick/IDS system. Com- 
parison with those galaxies will be used to assess the 
quality of our transformation to the Lick/IDS system. 
Because G93 has data for both central Reff/8 and and 
global {Reff/2) indices, we have a sample of 30 mea- 
surements. Our spectra were broadened to match the 
Lick/IDS resolution quoted by G93 (8.2 A FWHM), 
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Fig. 3. — The broadening corrections for the most relevant spec- 
tral indices compared to those of Trager et al. (1998), showing that 
our values are systematically lower. Figure (b) shows the compar- 
ison, with G93, of the velocity dispersions, where Acr = (Tours - 
fG93- Our values are systematically lower. 



while his data had corrections for both the broadening 
and the emission removed to recover the raw index val- 
ues. The index comparison shows that we have carried 
out a satisfactory transformation to the Lick/IDS sys- 
tem, with a negligible systematic discrepancy of 2% (< 
lowrs - Ig93 >), and a dispersion less than 8%. Figure 
2 shows this raw index comparison, with the dotted line 
representing equal values, and the average errors plotted. 

4.1.2. Error Determination and Broadening Correction 

The main source of error in the measurement of the 
line-strength indices is Poisson noise, which we calculate 
following the analytical approach of Cardiel et al. (1998). 
An almost irrelevant contribution comes from the error 
in the adopted zero point of the rotation curve, used to 
de-redshift the spectra. The redshift errors are estimated 
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Fk;. 4. — External comparison, with G93, of the most relevant 
spectral indices after the broadening and emission correction have 
been applied. The AIndex = Indexou^s - Indexggs. 



by measuring the radial velocities of each spectrum with 
respect to the three different rest frame templates: the 
stellar template and two SSP models corresponding to 
the RED and BLUE spectral intervals. Through this 
calculation, several sources of error are taken into ac- 
count, including template mismatch, spectral intervals, 
wavelength calibration, etc. Both Poisson and redshift 
errors are finally quadratically added to give the total 
errors of the indices. 

A broadening correction is applied to the raw spectral 
indices in order to remove the effect of the velocity dis- 
persion of the galaxy. This correction, a feature of the 
Lick approach, is calculated for each galaxy and spectral 
index. Our templates are the best matching SSP-Vaz99 
models, which can be customized for each galaxy by 
broadening them with both the instrumental and galaxy 
dispersions. For each spectral index, the broadening cor- 
rection is the ratio between the value of the index in the 
non-broadened versus the broadened template. Figure 
3(a) shows the comparison of our broadening corrections 
with those of Trager et al. (1998), where the largest 
discrepancy is found for the H/3 index. A similar dis- 
crepancy was also found in the Trager et al. (1998) vs. 
G93 comparison. This important multiplicative correc- 
tion amplifies any error in the velocity dispersion of the 
galaxies, introducing a new source of uncertainty in the 
determination of the spectral indices. We have taken this 
effect into account in the final error estimate, adding it 
quadratically to the errors of the indices and those of the 
emission correction. 

4.L3. External Comparison of the Spectral Indices 

Because this work utilizes two subsamplcs, a compari- 
son to external data is unnecessary. However, a compar- 
ison with the G93 results demonstrates that despite an 
initial good match of the uncorrected indices, the broad- 
ening correction is a source of systematic discrepancies. 
This is due to its capacity for amplifying the discrep- 
ancies existing in the ao determinations. As shown in 



Paper I, the value of <7q is rather sensitive to the choice 
of the spectral region used for its calculation. Figure 4 
shows the external comparison (our results vs. G93) of 
the relevant spectral indices, after the broadening and 
emission corrections. Our corrected spectral indices are 
systematically lower than those of G93, except for Mgb. 
For instance, H/3 is 10.6% smaller, while Fe5270 is 1.4%, 
with a slight systematic trend, and Fe5335 is 3.6% lower. 
Considering that the broadening correction increases the 
index value, the observed deficit can be interpreted as 
an under-correction of our indices. This is a consequence 
of the fact that both our broadening corrections and ctq 
values are systematically lower than those of G93. The 
striking discrepancy between the H/3 indices is due to 
the large disagreement between the H/3 correction curves. 
The external comparison, ours vs. G93, of the velocity 
dispersion is presented in Figure 3(b). On average, the 
velocity dispersion of our galaxies is 2 % smaller than 
those of G93, probably due to different spectral ranges 
used to calculate ctq, as discussed in Paper I. 

4.1.4. Age, Metallicity and [Mg/Fe] Determination with 

TMB03 

The stellar population parameters are extracted using 
the models developed by TMB03, which take the abun- 
dance ratio eff'ects into account. In order to break the 
age- metallicity degeneracy, we use the [MgFc]/ = (Mgb 
(0.72 X Fe5270 + 0.28 x Fe5335))i/2 index, defined by 
TMB03, and the H/3 index. The higher-order Balmer- 
line indices, H7 and H5, suffer less from emission con- 
tamination, but we have avoided their use because, as 
reported by Korn, Maraston & Thomas (2005), they arc 
not only more sensitive to metallicity than the H/3 in- 
dex, but also significantly more sensitive to the [Mg/Fe] 
ratio. The stellar population parameter extraction is car- 
ried out through an iterative process. In the first step, 
we derive the age and total metallicity from the [MgFe]/ 
vs. H/3 diagram with an [Mg/Fe] ratio arbitrarily set 
to 0.0. As shown by TMB03, both the [MgFe]/ and the 
H/3 indices have low sensitivity to the [Mg/Fe] ratio and, 
consequently, these first derived ages and metallicities 
are generally close to their convergence values. 

For the second step of the iteration, we construct a 
Mgb vs. (Fe) diagram, taking the previously derived age 
into account. From that diagram we extract the [Mg/Fe] 
ratio and proceed to the third step. We regenerate the 
[MgFe]/ vs. H/3 diagram, but for the new value of the 
[Mg/Fe] ratio and extract new ages and metallicities. 
The iterative process continues until the three extracted 
parameters converge to a stable value. The parameter 
determination is performed through linear interpolation 
in the points of the index-index grids. A considerable 
fraction (18.5 %) of the analyzed spectra show too low 
values of the H/3 indices which place them outside the 
grid. For those cases, extrapolation is needed to mea- 
sure the parameters, although it is worth noting that 
[Mg/Fe] varies slowly, with a typical rate of 0.002 dex 
per Gyr. An error of, say 5 Gyrs, would mean an in- 
crease/decrease of just 0.01 dex in [Mg/Fe]. Results for 
the three parameters and their errors, for the Reff/8 
aperture are shown in Table 2, columns 6, 7 and 8. Fig- 
ures 5(a) and (b) show the galaxy positions in two repre- 
sentative index-index grids. Figure 5(a) corresponds to 
an [Mg/Fe] = 0, while Figure 5(b) corresponds to the 
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median ago of 8 Gyrs. This representation is not en- 
tirely legitimate, because the correct practice would be 
to plot each galaxy in the grid corresponding to its spe- 
cific [Mg/Fe] and age. We note that this work employs 
a 2004 version of TMB03, which includes the response 
functions from Korn, Maraston & Thomas (2005) and 
supersedes in several aspects the original TMB03 based 
on Tripicco & BeU (1995). 

4.1.5. External Comparison of the Param.eters 

We have ten galaxies in common with Thomas et al. 
(2005), whose spectral indices were extracted from G93 
by these authors. The differences with their population 
parameters are a natural consequence of the above dis- 
cussed discrepancies between the spectral indices (see 
subsection 4.1.3). While their ages are, on average, 1.7 
Gyrs lower than ours, their mctallicities are larger by 0.14 
dex. For the [Mg/Fe] ratios, their larger (Fe) values and 
smaller age conspire to produce a negligible discrepancy 
with our values, well within the errors of the parameter. 
We have just three galaxies in common with Mendes de 
Oliveira et al. (2005), whose Fe5270, Fe5335 and Mgb 
indices are up to 4% larger than ours. Due to the use of 
the same stellar population analysis and TMB03 models, 
we expect the same pattern of parameter discrepancies 
found with Thomas et al. (2005). 

4.2. Age, Metallicity and [Mg/Fe] Determination: The 
MCGB05 Approach 

The MCGB05 approach is based on the grid of syn- 
thetic stellar spectra computed by Barbuy et al. (2003), 
which spans the range 4600-5600 A for stellar parameters 
4000 < Te// < 7000 K, 5.0 < log g < 0.0, -3.0 < 0.3, and 
two values of chemical abundance [a/Fe], 0.0 and -1-0.4 
dex (the a-elements considered are O, Mg, Si, S, Ar, Ca, 
Ne, Ti). MCGB05 present SSP models for the indices 
Mg2, Fe5270 and Fe5335, computed with improved co- 
efficients for the fitting functions given by Barbuy et al. 
(2003), and add a new set of spectra with [a/Fe] = -1-0.2. 
In the MCGB05 approach, the [Z/H] and [Mg/Fe] param- 
eters are calculated by interpolation (or extrapolation) of 
the Mg2 ys.{Fe) diagram corresponding to the age of the 
galaxy. Age values are taken from the prior TMB03 ap- 
proach, because no age-sensitive lines were computed in 
Barbuy et al. (2003). 

Although the TMB03 and MCGB05 approaches are 
apparently equivalent in some aspects, they show signifi- 
cant differences in their handling of the [a/Fe] abundance 
ratio. TMB03 corrects the empirical fitting functions of 
Worthey et al. (1994) with the response functions of 
Korn, Maraston & Thomas (2005) to account for the 
[a/Fe] enhancements, while the SSP models of MCGB05 
are based on synthetic stellar spectra whose [a/Fe] de- 
pendence is explicitly established along with the model 
computation. The results from this approach are pre- 
sented in Table 2 (columns 9 and 10). Figure 5(d) shows 
a typical metal- metal grid used for the [Z/H] and [Mg/Fe] 
calculations. It is worth mentioning that the superposi- 
tion of the whole sample serves only for qualitative pur- 
poses, as the grid is strictly valid only for those objects 
with ages of 8 Gyr. 

4.3. The Vaz99 Approach 



The main advantage of the Vaz99 approach is that the 
high resolution SSP models are adapted to the data, in- 
stead of the other way around. Therefore, there is no 
need for matching or broadening corrections. Both the 
cTinstrumentai and the cFgaUxy are taken into account by 
comparing the observed spectra with the a-broadened 
SSP models. As a drawback of this approach, non-solar 
[Mg/Fe] abundance ratios are not incorporated in the 
models. As with the previous approaches, it is equally 
important to carry out an accurate emission correction 
and, due to the different spectral resolution of the data, 
wc have repeated the double calculation of the correction 
for the H/3 index, following the procedures of Section 3. 

4.3.1. Age and Metallicity 

The Vaz99 SSP models have been constructed from 
the Jones' (1999) empirical stellar library, with high 
spectral resolution (1.8 A FWHM). They predict the 
Spectral Energy Distributions (SED) of a range of age 
and metallicity combinations, covering the intervals 0.7 
< log{Z/ZQ) < +0.2 and 1 to 17 Gyrs. These SEDs span 
just a BLUE (3856 - 4476 A ) and a RED ( 4795 - 5465 
A ) spectral interval. 

The Vaz99 stellar population analysis, which is 
markedly different from TMB03, proceeds as follows: 
first, the SSP model spectra are broadened from their 
original resolution (1.8 A FWHM) to the total resolu- 
tion of the spectrum of each particular galaxy. This 
atot resolution is the quadratic sum of our instrumental 
broadening aina ~ 133 km s^^ and the galactic veloc- 
ity dispersion agai taken from Table 2 (column 2). Once 
the broadened SSP models match the resolution of each 
particular galaxy, the comparison proceeds by measur- 
ing the spectral indices of the whole SSP model set. Our 
observed indices are then superimposed on the H/3 vs. 
[MgFe]' grid, constructed with the SSP values. The age 
and metallicity parameters are estimated by interpolat- 
ing among the SSP grid points. The SSP model grid 
is only valid for each agai, given that a-ms is a constant 
for the whole sample. Although it is not the case for 
the present study, it would be also possible to smooth 
a sample of galaxies to the same atot, in order to allow 
the use of a common grid. An extrapolation of the SSP 
model grid is sometimes needed to account for metallic- 
ities larger than the model upper limit [Z/H] = -1-0.2. 
Results for the central aperture Reff/8 are presented in 
Table 2, columns 3, 4 and 5. 

Errors for ages and mctallicities are also estimated with 
the grid. The error bars run from the locus of maximum 
to minimum index values. Although the error bars of 
the spectral indices are symmetric around their central 
values, this is not the case for the age-metallicity param- 
eters, because the SSP model grid is neither orthogonal 
nor equally spaced. 

4.3.2. Relative Abundances 

In the approach followed by Vaz99, one measures the 
abundance ratio of the galaxy relative to the scale-solar 
one, calculating the Mg and Fe abundances separately 
(see e.g. Vazdekis, Trujillo & Yamada 2004). Accord- 
ing to Tripicco & Bell (1995), the Mgb spectral index is 
predominantly sensitive to the a element Mg, while the 
combined index Fe3 is sensitive to Fe. Both the [Zms/H] 
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Fig. 5. — Tools from the three approaches used to determine the stellar population parameters in the present study. Grids (a) and (b), 
from the TMB03 models, allow the calculation of age, [Z/H] and [Mg/Fe] through an iterative process. Grid (a) is valid for [Mg/Fc]=0.0 
and grid (b) for 8.0 Gyrs. Panel (c) shows one of the grids used in the Vaz99 analysis, to derive the age and [Z/H]. This particular grid 
is only valid for galaxies with tjo = 170 km/sec. The MCGB05 analysis uses grid (a) to calculate the ages and grid (d) to measure the 
[Z/H] and [Mg/Fe]. The grid represented here is valid for 8.0 Gyrs. Although all the galaxies from the field (solid) and from HCGs (open) 
are shown, their positions in the parameter spax;e have no quantitative meaning, because these grids axe only valid for particular values of 
some parameters. The proper procedure, used in this study, requires the construction of customized grids for each galaxy. 



and [Zi?e/H] abundances are extracted from the H/3 vs. 
Mgb and H/3 vs. Fe3 diagrams and the [Mg/Fe] ratio 
is ([Zjvfg/H] - [Zi?e/H]). The Vaz99 simple approach is 
hampered by inaccuracies in the [ZMg/H] determination, 
which generally requires extrapolation of the model grid 
for values larger than the upper limit [Z/H] = + 0.2, 
especially for the most massive galaxies. Furthermore, 
the small departures from orthogonality in the H/3 vs. 
metal-index grids transform errors in age into errors in 
metallicity. They have been incorporated into the errors 
obtained for the [Mg/Fe] ratios. 



4.4. Comparison of the Three Approaches 

It is worth making a comparison of the extracted pa- 
rameters in order to evaluate the systematic effects in- 
troduced by the choice of the model and stellar pop- 
ulation analysis. This is also interesting, outside the 
context of the present study, in order to better inter- 
pret other stellar population results obtained using any 
of these different approaches. In the present study, the 
TMB03 approach has been chosen as the common refer- 
ence for comparison of all our results, for the sole reason 
that it makes comparison with the literature easier, as 
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explained in section 5.2. Figure 6 shows the comparison 
with Vaz99, while Figure 7 shows the comparison with 
MCGB05. 

4.4.1. The Vaz99 vs. TMB03 results 

As shown in Figures 6, ages obtained with the TMB03 
analysis are, on average, 16% lower than those obtained 
with the Vaz99 approach, while TMB03 metallicitics arc 
0.16 dex larger, on average, than those of Vaz99. These 
discrepancies can be qualitatively explained by differ- 
ences between the TMB03 and Vaz99 model grids, as 
shown in Figure 6(b). In order to make the comparison 
easier, all the grids of Figure 6 have been adapted to 
show the same age and metallicity points. Note that it is 
inaccurate to represent both TMB03 and Vaz99 grids in 
a common graph, except for the special case where the 
observed crtot = {'^lai + '^InstY^'^ equals the instrumen- 
tal broadening of auck = 200 km s~^. For our particular 
(Jinsti this means that both grids are roughly equivalent 
when (jgai = 172 km s~^. In Figure 6(b), we have used 
the Vaz99 grid for a agai = 172 km s^^ galaxy and the 
TMB03 grid with [Mg/Fe] = 0.2, the average value of 
the sample. From Figure 6(b) we can understand qual- 
itatively the observed age and metallicity discrepancies. 
For instance, a galaxy (solid square) with agai = 172 
km , H/3 = 1.72 and [MgFc]/= 3.48 would have an 
Ageyaz99 = 11-0 Gyrs and an AgCTMBoa = 9.0 Gyrs, as 
well as [Z/R]vaz99 = 0.0 and [Z/HJtmbos = 0.17. The 
relatively lower dispersion of the points in Figure 6(c) 
with respect to (a) is probably due to the fact that lines 
with cqual-metallicity are more parallel than lines with 
equal-age. 

Concerning the abundance ratios, the [Mg/Fe] ob- 
tained with the Vaz99 approach are larger than the 
[Mg/Fe] of TMB03. Furthermore, the discrepancy in- 
creases with the overabundance. These differences might 
be attributed to the fact that, according to Korn, Maras- 
ton & Thomas (2005), Mg and Fe indices show inverse 
sensitivities to [Mg/Fe], i.e. the Mgb index increases 
with [Mg/Fe], while the Fe3 index decreases. As these re- 
sponses are used in TMB03 models, the Vaz99 approach, 
which is scaled-solar, consequently overestimates [Mg/H] 
and underestimates [Fe/H] metallicitics with respect to 
TMB03, leading to systematically larger [Mg/Fe] values. 
A similar graphic explanation of the discrepancy is pre- 
sented in Figures 6(d) and (f), which show the H/3 vs. 
Mgb and H/3 vs. Fe3 grids for [Mg/Fe] = 0.0 (Vaz99) 
and [Mg/Fe] = 0.4 (TMB03). In the Vaz99 approach, 
both the [Mg/H] and the [Fe/H] metallicitics are mea- 
sured with respect to the [Mg/Fe] = 0.0 grid (solid-line). 
For any [Mg/Fe] > 0, the Vaz99 approach overestimates 
[Mg/H] (Figure 6(d)) and underestimates [Fe/H] (Fig- 
ure 6(f)), compared to TMB03, and leads to the ob- 
served larger [Mg/Fe] (Vaz99) values. In order to quantify 
the discrepancy, we have represented an arbitrary galaxy 
(solid triangle) with agai = 172 km in Figures 6(d) 
and (f). According to TMB03, this galaxy would have 
[Mg/Fe] = 0.4, because both [Fe/H] and [Mg/H] metal- 
licitics converge to the same value (0.0) only for this par- 
ticular [Mg/Fe] ratio. In the Vaz99 approach, the same 
galaxy would have [Mg/Fe] = 0.76, because [Mg/H] = 
0.26 (Figure 6(d)) and [Fe/H]= - 0.50 (Figure 6(f)). This 
is the typical disagreement observed in Figure 6(e) for the 



objects with the largest [Mg/Fe] ratios. 

4.4.2. The MCGB05 vs. TMB03 results 

Comparison of MCGB05 with TMB03 results, pre- 
sented in Figures 7 (a) and (c) shows that the MCGB05 
metallicitics are, on average, 0.07 dex lower than TMB03. 
The [Mg/Fe] ratio is systematically larger in MCGB05, 
with the discrepancy increasing with the abundance ra- 
tio. The comparison is now easier to interpret because 
both methods share the same age (TMB03) values. The 
Mg vs. Fe diagrams for each approach are presented in 
Figures 7(b) and 7(d), where a common age (10 Gyr) 
has been selected. Two arbitrary galaxies have been 
represented and their [Z/H] and [Mg/Fe] parameters 
have been extracted through inter- or extrapolation of 
the grid values. In order to guarantee the equivalence 
of the galaxies, we have matched the two Mg indices 
with the relation Mg;, = 12.4 x Mg2 + 0.91. This is 
a low-dispersion relation obtained empirically from our 
data, including 80 different measurements. The first ex- 
ample, the empty triangle, with [Z/H]tmbo3 = +0.05 
and [Mg/Fe]TMB03 = 0.10 shows comparable parame- 
ter values in MCGB05, with a slightly lower [Z/H] (Fig- 
ure 7(d)), while the second one, the solid circle, with 
[Z/H]tmbo3 = +0.23 and [Mg/Fe] TMBos — 0.40 shows 
an inflated abundance ratio (0.66), as expected. 

In order to provide a common result for the three ap- 
proaches, the Vaz99 and the MCGB05 partial results 
have been transformed into the TMB03 common frame 
of reference, by correcting their mutual systematic dis- 
crepancies. From the linear fits obtained from Figures 6 
and 7 we get: 

AgeTMBoa = 0.87 x Agevaz99 - 0.27 
[Z/H]tmbo3 = 0-75 X [Z/H]y^^,, + 0.17 
[Mg/Fe] 

TMB03 ~ 

0.33 X [Mff/Fe]^„,99 + 0.10 

[^/^Itmbos =0.56 X [Z/H]^cGB05 + 0-07 
[Mg/Fe]^j^so3 = OA8 X [Mg/Fe]j^cGB05 + 0-07 



The final results of the present study, included in Table 
2 (columns 11, 12 and 13), are an average of the individ- 
ual results from the three approaches, after they have 
been converted to the common TMB03 frame. The fi- 
nal errors are computed from the errors of the individual 
approaches added in quadrature. 

From Figures 6(c) (e) and 7(a) (c) it is easy to con- 
clude that the metallicity parameters from the Vaz99 
and MCGB05 models are in better mutual agreement 
than those of TMB03. This provides a sound consistency 
test for these two widely different approaches (Vaz99 and 
MCGB05). 

5. RESULTS 

We expect the ages, metallicitics and [Mg/Fc] abun- 
dance ratios of the galaxies to vary with both the inter- 
nal and environmental characteristics. We have consid- 
ered the galaxy mass, or its related observable the central 
velocity dispersion ctq, as the most appropriate internal 
parameter. Concerning the environment, although our 
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Fig. 6. — Comparison of the stellar population parameters obtained with the TMB03 and Vaz99 models. Panels (a), (c) and (e) show 
the comparison of ages, [Z/H] and [Mg/Fe] for the whole sample, where systematic discrepancies are evident. Panels on the right provide 
qualitative explanations for the discrepancies. In figure (b), equivalent grids from both models have been superimposed. For the same pair 
of spectral indices, an arbitrary galaxy (the solid square) shows systematically larger age and lower [Z/H] with the Vaz99-grid (continuous 
line) than with the TMB03-grid (dashed line), in agreement with panels (a) and (c). The systematic discrepancy in [Mg/Fc] is explained 
in figures (d) and (f) for an arbitrary galaxy (solid triangle) with [Mg/Fe] (TMB03) = 0.4, i.e. a galaxy with the same [Mg/H]=[Fe/H]=0.0 
in the [Mg/Fe]=0.4 grids of the TMB03 models (dashed lines). The Vaz99 procedure to measure [Mg/Fe], i.e. [Mg/Fe] =[Mg/H]-[Fe/H], 
gives a much larger value ([Mg/Fe] = 0.26 - (-0.50) = 0.76), as observed in panel (e). 
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Fig. 7. — Comparison of stellar population parameters obtained with TMB03 and MCGB05 models. Left panels show the comparison 
between [Z/H] and [Mg/Fe] with the equal- value line (dot-dashed) and the least-squares fit (dotted) to show the systematic discrepancies. 
Panels (b) and (d) explain qualitatively the discrepancies, by showing two arbitrary galaxies with [Mg/Fe] (TMB03) = 0.1 (open triangle) 
and with [Mg/Fe] (TMB03) = 0.4 (solid circle). These same galaxies, when studied with the MCGB05 models, have [Mg/Fe] (MCGB05) = 
0.1 and 0.66, respectively, as observed in (c). The [Z/H] shows a smaller discrepancy, as seen in (a). 
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study is based in the simplistic division into HCG and 
the field, wc have also explored the diversity of the het- 
erogeneous HCG environment by means of the crossing 
time which characterize their dynamical evolution. Ob- 
servations of the same galaxy with different slit inclina- 
tions are not treated as different objects; instead only 
the observation with the highest S/N is considered. 
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5.1. Stellar Population Parameters vs. Velocity 
Dispersion 

The observable parameter ctq can be indirectly linked 
to the stellar mass through the mass-to-light ratios and 
the Faber- Jackson relation (e.g. equation (2) in Thomas 
et al, 2005). In Figure 8 we show the correlation of the 
stellar population parameters with ctq for the two dif- 
ferent environments. Data from the literature have also 
been included in the graphs and used to calculate the 
Icast-squarcs fits, as described in subsection 5.2. The 
dwarf galaxy HCG 37e has been excluded from the fits 
of Figures 8(c) and 8(f) to avoid its excessive statistical 
weight. Although only our average results are presented 
in Figure 8, the same qualitative behavior has been con- 
sistently found using the results from each of our three 
alternative approaches (sec Table 2). 

Several of the trends found in Figure 8 deserve com- 
ment. Although our data do not show significant cor- 
relation between age and velocity dispersion (the solid 
symbols in Figure 8(a) and (b)), they follow the large- 
dispersion trend for field data noted in the literature 
(Trager et al. 2000a; Thomas et al. 2005), and in agree- 
ment with the extensive work of Gallazzi et al. (2006) 
using SDSS early-type galaxies. Both the fit to all the 
data (dotted line) and the median value of only our re- 
sults (dashed line) are shown in Figure 8(b). For galaxies 
in HCGs (Figure 8(a)), ages do not show any trend with 
log Co and only the median value of our results is plotted. 
There is a positive correlation between the [Z/H] metal- 
licity and log (Tq for the galaxies in both environments 
(Figure 8(c) and (d)), which is remarkably steeper in the 
groups. A positive correlation is generally interpreted 
as the action of deeper potential wells hindering the loss 
of metal-enriched gas through galactic winds (Greggio 
1997). We interpret the steeper trend of the [Z/H] vs. 
log ctq in the group environment as the action of a SF 
truncation mechanism which, nevertheless, does not pre- 
clude the aforementioned processes. 

The field galaxies show a clear correlation between 
[Mg/Fe] and log (Tq (Figure 8(f)), in good agreement with 
the data from the literature. This trend has been inter- 
preted as a sign of downsizing, meaning that the more 
massive galaxies complete their star formation faster (e.g. 
Thomas et al 2005). A novel result of the present study 
is the negative correlation found between [Mg/Fe] and log 
(To of the low-mass galaxies of HCGs in contrast with the 
behavior of the corresponding field galaxies. Figure 8(e) 
shows the two contrasting linear fits, with a dashed line 
for the HCGs and a dotted one for the field. Despite the 
striking visual contrast between both trends, the size of 
the [Mg/Fe] errors make the difference between slopes 
barely significant, just slightly above the onc-o" level. 
Nevertheless, other peculiar behaviors reinforce the re- 
ality of the discrepancy, like the deficient metallicity of 
low-mass galaxies in Figure 8(c) and other photomet- 
ric and environmental parameters discussed below. We 
have arbitrarily selected a subsample, including the low- 
mass HCG galaxies with 2.0 < logao < 2.2, because they 
are largely responsible for the discrepancy between the 
slopes of the fits in Figure 8(e). At the same time, mem- 
bers of this low-mass subsample, marked with circles in 
Figure 8, also populate a region of the parameter space 
where [Mg/Fe] > 0.21 and [Z/H] < 0.16. Hereafter, we 



call them anomalous galaxies, noting that this anomaly 
refers only to their discrepancy with the field galaxies. 
Dwarf galaxy HCG37e is excluded from the subsample, 
as its behavior coincides with that of the field galaxies. 
In agreement with the downsizing concept, this galaxy 
presents parameter values comparable to the Sun, sug- 
gesting an extended SF history, similar to our galaxy 
disk. 

There are some remaining sources of concerns regard- 
ing the reality of the anomalous enhancement of [Mg/Fe] 
in the low-mass HCG galaxies. Among them, first is the 
possible bias introduced by the emission correction. As 
explained in subsection 5.4.1, the anomalous galaxies are 
preferentially AGNs, i.e. they are contaminated by emis- 
sion. Errors in the emission correction are propagated 
unevenly to the stellar population parameters. Their 
sensitivity to changes in H/3 is largest for the age and 
lowest for the [Mg/Fe]. For instance, an unrealistically 
large error in the correction of AH/3 ~ 0.3 is converted 
into Aagc - 7 Gyrs, A[Z/H] ~ 0.2 dex and A[Mg/Fe] 
~ 0.02. Figures 5(a) or 5(c) are useful to understand the 
qualitative effect of a change in H/3. Only a scenario in 
which the emission-contaminated spectra of the low-mass 
galaxies are systematically under-corrected by an exag- 
gerated amount would reproduce the behavior observed 
in the [Z/H] vs. log ctq trend of the HCG galaxies, al- 
though it would leave the anomalous trend in [Mg/Fe] 
vs. log ctq almost unaffected. Considering the double 
emission correction carried out in section 3, we have no 
reason to suspect such a tweaked under-correction. 

A further source of concern is that this anomalous ef- 
fect has passed unnoticed in previous studies, i.e. Proctor 
et al. (2004) and Mendes de Oliveira et al. (2005). As 
we show in the next subsection, their results generally 
agree with ours, but the shortage of data for low-mass 
galaxies, in the interval logao = 2.0 - 2.2, prevented these 
authors from observing the effect. 

5.2. Data from the Literature 

Several other studies on the stellar populations of 
early-type galaxies can be used to assess the correlations 
found here and to expand our samples. For field galax- 
ies, we have used the works of Thomas et al. (2005) and 
Trager et al. (2000a). Least-squares fits to the correla- 
tions of the stellar population parameters, i.e. Age, [Z/H] 
and [Mg/Fe], with logcro have been calculated for our 
data and the external samples. Despite slight systematic 
discrepancies, explained below, the average agreement 
among the three samples is always within the rms devi- 
ations around the fits, demonstrating good agreement of 
our field galaxy results with the literature. 

In addition to providing indpendent comparisons, the 
works of Thomas et al. (2005) and Trager et al. (2000a) 
can be used to expand our sample, because they share 
several galaxies in common with our field subsample, 
which allow for conversion to the TMB03 reference frame 
of our results. Galaxies from the "low-density" subsam- 
ple of Thomas et al. (2005) have been corrected to our 
reference frame; via 12 common galaxicis. The conver- 
sions for the stellar-population parameters, a measure 
of the systematic discrepancies, are Agecorr = 1-07 x 
AgCT/iomas + 1-24; [Z/H] 

corr — 

0.69 X [Z/H] 

Thomas 

0.02 and [Mg/Fe]corr = 1-1 x [Mg/Fe]Thomas- After ex- 
cluding the galaxies which are redundant with the Trager 
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Fig. 8. — Correlations of our results with logcro, comparing the behavior of elliptical galaxies in the FIELD (right-hand panels) with their 
counterparts in Hickson Compact Groups (HCG) (left-hand panels). Data from the literature (open symbols) have been superimposed to 
assess our results. They were taken from Thomas et al. (2005, T05), Trager et al. (2000a, TOO) and Mendes de Oliveira et al. (2005, 
MCGB05) and previously converted to our reference frame. Least-squares fits (dotted lines) for all the field galaxies are displayed in figures 
(b), (d) and (f) and repeated in the group galaxy panels (c) and (e). Dashed lines in figures (a) and (b) show the median value of our 
data, while panels (c) and (e) show the least-squares fit (dashed) to all the objects, excluding the galaxy HCG37e, with log(To=1.93. The 
anomalous galaxies (see text) are tagged with a circle and the average error bars of our data are plotted in the low-right corner of each 
panel. 



et al. (2000a) sample, 41 field galaxies are plotted in 
Figures 8(b), (d) and (f). A similar procedure has been 
used to convert the Trager ct al. (2000a) parameters to 
our reference frame, using the relations AgCcorr = 0.95 
X A 1.76; [Z/R]corr = 0.27 x [Z/H] 

1 rager i 

0.15 and [Mg/Fe]corr = 1-05 X [Mg/Fe]Trager + 0.05. 



The resulting sample, comprising 40 field galaxies, is also 
presented in Figures 8(b), (d) and (f) where a common 
least-squares fits to all the data, including ours, have 
been displayed. 

Proctor et al. (2004) and Mendes de Oliveira et 
al. (2005) studied elliptical galaxies in compact groups 



16 



which can be used to assess our resuhs. The first study 
is especially difficult to convert to our reference frame, 
because it follows a unique approach, very dissimilar to 
ours, and has only three galaxies in common. In ad- 
dition, we suspect that their [Mg/Fe] results differ sys- 
tematically from ours because their field galaxies do not 
show the striking [Mg/Fc] vs. logcro correlation seen in 
our data and those of other authors, as displayed in Fig- 
ure 8(f)). The comparison with Mcndes de Oliveira et 
al. (2005) is simpler because they use the same TMB03 
stellar population analysis. The conversion is carried out 
in two steps, starting with a transformation of their data 
from an aperture of Re///2 to Re///8, using the average 
parameter gradients obtained from our data, erg = 1-1 x 
cr2 - 14.8. Because the reference frame of the Mendes 
de Oliveira et al. (2005) data is the G93 sample reduced 
with the TMB03 models, we use our twelve galaxies from 
the G93 sample to derive the conversions. The final rela- 
tions are Ucorr = 0.96 X UMendes + 24.3; kgCcorr = 1-04 
X AgeMendes + 0.82; [Z/H]corr = 0.49 X [Z/H]Me«des + 

0.04 and [Mg/Fe]corr = 0.94 x [Mg/¥e]Mendes + 0.05, 
with the results displayed in Figures 8(a), (c) and (e). 

The conversion of the Mendes de Oliveira et al. (2005) 
data to our TMB03 frame adds a new object to the crit- 
ical interval logcro = 2.0 - 2.2 where the anomaly is de- 
tected. Their galaxy HCG76d, at logcro = 2.19, shows 
the typical anomalous behavior, with enhanced [Mg/Fe] 
= 0.47 and reduced [Z/H] = -f0.07. 

5.3. Stellar Population Parameters vs. Environment 

Prom our simplistic environmental division into HCGs 

and the field, we can conclude from Figure 8 that our 
galaxies in the HCGs are older and more metal-poor than 
their counterparts in the field, in agreement with previ- 
ous works, e.g. de la Rosa et al. (2001b), Proctor et al. 
(2004) and Mendes de Oliveira et al. (2005). Our galax- 
ies in groups are, on average, 1.6 Gyrs older than in the 
field. When field galaxies are discriminated according 
to their local density, those located in the less dense en- 
vironments are the youngest. For instance, of the five 
oldest galaxies in Figure 8(b), four are located in high 
density environments, according to the classification of 
Thomas et al. (2005). Galaxies NGC4552, NGC4649 
and NGC4697 are located in the Virgo cluster, while 
NGG7619 belongs to the Pegasus-I group. Would they 
have been excluded from the sample, the difference be- 
tween groups and the field (low density) would increase 
to 4.3 Gyrs. The metallicity [Z/H] values for the field 
galaxies are 0.11 dex larger, on average, than those of 
the CGs, while the average [Mg/Fe] ratios are identical 
in both environments. Much more interesting environ- 
mental dependencies of the [Z/H] and [Mg/Fc] have been 
mentioned in the previous subsection, concerning their 
trends with log ctq. 

From a dynamical point of view, HCGs form a hetero- 
geneous sample, containing groups at different evolution- 
ary stages (Ribeiro et al. 1998). A convenient parameter 
to characterize the degree of dynamical evolution of the 
group is the crossing time (Hotc), whose values we have 
taken from Hickson et al. (1992). For instance, these au- 
thors found a deficit of late-type galaxies in groups with 
shorter tc, where galaxy-galaxy interactions are supposed 
to be more frequent. 

In our study, Figures 9(a) and (b), we can see that 
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Fig. 9. — Dependence of the parameters for galaxies in HCGs 
with the dynamical state of the group, represented by the crossing 
time log(Hotc). Least-squares fits (dashed lines) are displayed in 
panels (a) and (b), with the average value (dot-dashed lines) of the 
corresponding parameter in the field, except for panel (c), where 
no correlation has been found for [Mg/Fe]. Anomalous galaxies, 
defined in the text, have been identified with a circle. 



both age and metallicity depend on the dynamical evo- 
lution of the HCGs. Least-squares fits arc drawn to show 
that stellar population parameters for the less dynami- 
cally evolved HCGs, i.e. high log(Hotc), closely resemble 
those of the field, while the parameters of the HCGs with 
a small log(Hotc) show the largest disagreement with the 
field. This is a reasonable result, as it is generally ac- 
cepted that galaxy interactions are the main driver of 
the environmental differences. The [Mg/Fe] does not 
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TABLE 3 

Peculiarities of the anomalous galaxies 



Subsample ao [Z/H] [Mg/Fe] log(Hotc) Emission a4>0 



Anomalous 146.6 ± 18.4 -0.02 ± 0.15 0.30 ± 0.08 -2.01 ± 0.19 67 % 100 % 

Rest 207.4 ± 24.5 0.15 ± 0.16 0.21 ± 0.06 -1.82 ± 0.22 33 % 46 % 



show a clear dependence on the crossing time, (Figure 
9(c)), but it is worth noticing that the anomalous galax- 
ies, marked with circles, are preferentially located in the 
more evolved HCGs (log(Hotc) < -2.0, suggesting that 
the source of the anomaly is related to the frequency of 
galaxy interactions. 

5.4. Which are the "Anomalous" Galaxies? 

The anomalous galaxies in our sample are HCG14b, 
HCG46C, HCG59b, HCGlQa, HCG46a and HCG57f, 

although the first three are the most extreme cases. 
From external samples, we can add to this list HCG76d 
(Mendes de Oliveira et al. 2005). Table 3 shows a com- 
parison of the properties of the anomalous with the re- 
mainder of the HCG galaxies in our sample, excluding 
the dwarf HCG37e. As already mentioned, the anoma- 
lous galaxies show, on average, lower ctq (col. 2) and 
[Z/H] (col.3) as well as a larger [Mg/Fe] (col. 4) than 
the rest of the HCG sample. They are also located in 
more dynamically evolved groups, having shorter cross- 
ing times (col. 5). We can estimate the stellar masses of 
the galaxies based on the relation given by Thomas et al. 
(2005), logM, « 0.63 + 4.52 loga. The average mass of 
the anomalous galaxies thus calculated is 2.6xlO^°M0, 
while for the rest of the sample it is I.SxIO^^Mq. 

5.4.1. Activity 

Two thirds of the anomalous galaxies show nebular 
emission, against only one third of the rest of the HCG 
galaxies. We have studied their ionization mechanisms 
with the standard procedure, based on the characteristic 
spectral line ratios of the different activity classes (e.g. 
Coziol et al 1998). The emission lines were measured in 
the template-subtracted spectra, by following the same 
procedure described in section 3. According to the stan- 
dard diagnostic diagram (Veilleux & Osterbrock 1987), 
which uses the [Nil] /Ha and [0III]/H/3 line ratios, only 
one of the emission-line galaxies shows weak star for- 
mation activity, while the rest are classified as AGNs. 
All active galaxies increase their [OIII]A5007 emission in- 
wards, from Reff/2 to Reff/8, reinforcing the idea of a 
nuclear source of activity. 

5.4.2. Photometric Profiles 

Mendes de Oliveira &: Hickson (1994) find that three 

of our anomalous galaxies, namely HCG 14b, 57f and 
59b, stand out for having extended profiles and being 
more diffuse than other galaxies of the same luminos- 
ity. They note that these galaxies are similar to the cD 
galaxies in clusters, albeit much less luminous. Isophotal 
photometry for the galaxies in this sample (Ferreras et al. 
2006) shows that all the anomalous galaxies with mea- 
sured a4 parameters have disky isophotes (positive a4), 



unlike the rest of the sample, where 54% of the galaxies 
have boxy profiles (negative a4). This can be a selec- 
tion effect, as disky galaxies tend to be less luminous 
than their boxy counterparts, although with a substan- 
tial amount of crossover (e.g. Bender et al. 1989). 

6. DISCUSSION 

Unlike the largest ellipticals, which are already inactive 
by z « 1, medium-sized and dwarf galaxies in the field are 
still actively forming stars in the present universe. This 
phenomenon, known as "downsizing", has been repeat- 
edly observed (Cowie et al. 1996; Kauffmann et al. 2003; 
Thomas et al, 2005) and recently interpreted in terms of a 
feedback-regulated picture of galaxy formation (Scanna- 
pieco. Silk & Bouwens, 2005). Studies based on redshift 
surveys comprising vast samples of galaxies show a clear 
bimodality in star formation properties, with two dis- 
tinct galaxy populations (Baldry et al. 2004). The red 
sequence comprises passive galaxies without SF, while 
the blue sequence includes galaxies with active SF. The 
global decrease observed in the SF rate density since red- 
shift 1-2 (e.g. Hopkins & Beacom 2006) suggests that 
galaxies are systematically transferred from the blue to 
the red sequence. Several observations can be used to 
outline the properties of the mechanism responsible for 
the transition or SF-truncation. The mere existence of 
the bimodality implies that the mechanism must be fast 
in order to preserve the separation of the sequences. A 
slow SF truncation mechanism is not viable, as it would 
smear out that separation. According to a study by Ma- 
teus et al. (2006), the galaxies filling the gap between 
both color sequences are mainly AGN hosts, indicating 
that the SF truncation mechanism must also involve nu- 
clear activity. Concerning the preferred environment in 
which the mechanism might act, it has been repeatedly 
found that dense environments are hostile to SF (e.g. 
Baldry et al. 2004). However, extreme densities are not 
needed, because Lewis et al. (2002) and Gomez et al. 
(2003) found that the SF truncation mechanism operates 
at the low local densities found well beyond the virial- 
ized cluster cores. Mateus & Sodrc (2004) went further, 
showing that the SF in galaxies is affected by a wide 
range of environments, including galaxy groups and the 
field. A third characteristic of the SF truncation mech- 
anism is that it must be able to act everywhere. Finally, 
Kauffmann et al. (2003) found that the transition be- 
tween both sequences happens at a critical stellar mass 
3x1O^"M0, suggesting that this is the typical mass of 
the galaxies undergoing SF-truncation. 

Although there is no shortage of mechanisms capable 
of truncating star formation in galaxies, none of them ful- 
fills the three expected characteristics. For instance, ram 
pressure stripping (Dressier & Gunn 1983) needs high ve- 
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locities and densities at the inter-galactic medium, which 
are only found in dense clusters. The perturbation of the 
gas-rich halo of the galaxy, called strangulation (Larson, 
Tinsley & Caldwell 1980) produces a too slow decrease 
of the SF, in > 2 Gyrs (Balogh et al. 2004). Never- 
theless, recent hydrodynamical simulations (Di Matteo, 
Springel & Hernquist 2005) have shown that galaxy m,erg- 
ers, commonly considered a mechanism for SF enhance- 
ment, possesses powerful SF truncation properties. Ac- 
cording to the simulations, the galaxy merger induces 
a short-lived star formation burst and the triggering of 
an AGN, which rapidly depletes the gas and quenches 
the star formation. Consequently, the merger of galaxies 
fulfills the three characteristics required for the SF trun- 
cation mechanism, i.e. it is fast, involves AGNs and can 
potentially act everywhere. 

Compact groiips of galaxies, with their high densities 
and low velocities make the ideal environment for merg- 
ers to occur. This line of reasoning lead us to concen- 
trate our search for SF-truncation tracers in this minority 
class. Considering the plethora of faint galaxies detected 
in HCGs (de Carvalho et al. 1997), it is worth consid- 
ering the possibility of minor mergers involving satellite 
or dwarf companions (Hernquist & Mihos 1995; Bell et 
al. 2005). Galaxy-galaxy interactions, and especially 
minor mergers, are expected to occur in loose groups or 
even the field. In summary, we propose a scenario in 
which galaxies that enter the group environment suffer a 
minor merger. The consequences of the merger depend 
on the mass of the entering galaxy. For instance, due 
to the downsizing, massive galaxies have completed their 
SF long before entering the group and therefore only suf- 
fer mild consequences. A dwarf galaxy, with an extended 
SF history, will just contaminate the massive galaxy with 
[Z/H]-ricli and [Mg/Fe]-poor stellar populations. How- 
ever, a massive galaxy is passive and there is no SF to 
truncate. In contrast, a minor merger with a low-mass 
galaxy would likely induce more drastic consequences. 
To begin with, an unperturbed low-mass galaxy is ex- 
pected to be still actively forming stars in the present 
universe. The merger event would deplete its gas and 
truncate any further SF, leaving the values of the stellar 
population parameters frozen. Unperturbed in the field, 
its counterparts with similar mass steadily increase their 
metallicity [Z/H] and decrease the [Mg/Fe], by incorpo- 
rating Fe in the newly formed stars. 

Our observations fit this scenario well. First of all, 
we observe that galaxies in CGs are, on average, older 
and more metal-poor than their counterparts in the field, 
reminding us of the hostility to SF which characterizes 
dense environments. The anomalous enhanced [Mg/Fe] 
abundance ratio and deficient metallicity detected in 
the low-mass galaxies of our HCG sample, can be in- 
terpreted in the frame of the aforementioned merger 
+ SF-truncation mechanism. We emphasize that both 
the [Mg/Fe] ratio and the [Z/H] are not constant but 
evolve along the SF history of the galaxy. Soon after 
the SF history starts, the SN-II enrich the ISM with a 
elements which are readily incorporated into the newly 
forming stars. At that moment, the [Mg/Fe] ratio at- 
tains its highest value, which later decreases as soon as 
the Fe is released by the SNIa and incorporated into the 
stars. Similarly, along the SF history, the global metal- 
licity [Z/H] of the stellar population increases. The SF- 



truncation puts an end to this evolution and leaves the 
galaxy with their parameter values at the time of the 
merger event. 

The anomalous galaxies with larger [Mg/Fe] and lower 
[Z/H], are simply galaxies which truncated their SF. The 
fact that they belong to the densest groups (shorter 

cross), show larger incidence of AGN and a predomi- 
nance of disky profiles reinforces the idea that their SF 
was truncated by a minor merger. Simulations show that 
disky ellipticals are the result of mergers between pro- 
genitors of disparate masses (Khochfar & Burkert 2005; 
Jesseit, Naab & Burkert 2005). The fact that our anoma- 
lous galaxies have an average stellar mass of 2.6x 10^°Mq 
locates them in the transition region between the red 
and blue sequences defined by Kauffmann et al. (2003), 
where the SF-truncation is likely taking place. 

7. SUMMARY AND CONCLUSIONS 

In this paper we have studied the stellar populations of 
elliptical galaxies in HCGs in order to compare them with 
their counterparts in the field. Assuming that galaxies in 
groups have been incorporated from the field, differences 
between members of those environments would unmask 
the long sought mechanism which makes the dense envi- 
ronments hostile to star formation. 

Our stellar population study has tried to minimize 
the sources of systematic effects by carrying out several 
parallel stellar population analyses with different pro- 
cedures, models and emission corrections. Although the 
three alternative models provide different absolute values 
for the stellar population parameters, i.e. age, metallicity 
and abundance ratio, their results, concerning the differ- 
ences between both field and HCG galaxies are qualita- 
tively similar. 

Prom our results we can extract the following conclu- 
sions: 

(i) Galaxies in the HCGs are, on average, 1.6 Gyrs 
older and more metal-poor, in 0.11 dex of [Z/H], than 
their counterparts in the field. This result is in agreement 
with the previous findings reported in the literature and 
supports the view of a hostility of the dense environments 
to star formation. 

(ii) The [Mg/Fe] of galaxies in the field show a trend 
with their velocity dispersion (ctq). More massive galax- 
ies (larger ctq) show larger abundance ratios. This result, 
also in agreement with previous studies, has been inter- 
preted as a sign of downsizing, where the less massive 
galaxies have more extended SF histories. 

(iii) The low-mass galaxies in HCGs (with average 
M^. « 2.6xlO^°M0) show a distinct behavior with re- 
spect to similar galaxies in the field. Those in HCGs 
show higher [Mg/Fe] abundance ratios and lower [Z/H] 
metallicities. These anomalous galaxies also show other 
peculiarities when compared to the rest of the HCG sam- 
ple. They are preferentially found in the most dynami- 
cally evolved HCGs (lower crossing times) and show an 
excess of AGNs, disky isophotes and extended photomet- 
ric profiles, reminiscent of those in cD galaxies. 

We propose that low-mass galaxies from the field, ac- 
tively forming stars, are incorporated into the dense CG 
environment, where they suffer mergers with other galax- 
ies. As a result of the merger, the AGN feedback depletes 
the gas and quenches star formation. This SF truncation 
blocks the natural decrease of the [Mg/Fe] which accom- 
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panies any extended star formation. It equally hinders 
the metal enrichment and leaves a galaxy with larger 
[Mg/Fe] and lower [Z/H], as observed here. Massive 
galaxies, on the contrary, are excluded from this behav- 
ior because they are passive in the present universe and 
there is no SF to truncate. In summary, we conclude that 
we have detected traces of the merger action in low-mass 
galaxies of HCGs. 

If the merger + SF-truncation mechanism is accepted 
to be working in the Compact Groups, it can contribute 
to understanding a long standing apparent contradiction 
in HCGs, where their galaxies, despite showing abundant 
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signs of tidal interactions, do not show the expected star 
formation enhancement. 

The merger + SF-truncation would possibly be the 
mechanism responsible for the global decay of the SF rate 
since redshift 1-2, as the groups are the places where the 
galaxies are pre-processed and transformed into passive 
objects, before being incorporated into clusters. 
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